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(^ This paper presents a new formulation for trailing edge noise radiation from 

m rotating blades based on an analytical solution of the convective wave equation. 

^_^ It accounts for distributed loading and the effect of mean flow and spanwise 

^ wavenumber. A commonly used theory due to Schlinker and Amiet (1981) predicts 

J>~> trailing edge noise radiation from rotating blades. However, different versions of 

"^ the theory exist; it is not known which version is the correct one, and what 

^ the range of validity of the theory is. This paper addresses both questions by 

rri deriving Schlinker and Amiet 's theory in a simple way and by comparing it to the 

^ new formulation, using model blade elements representative of a wind turbine, a 

Q cooling fan and an aircraft propeller. The correct form of Schlinker and Amiet's 

■^ theory (1981) is identified. It is valid at high enough frequency, i.e. for a Helmholtz 

^»-> number relative to chord greater than one and a rotational frequency much smaller 

i-G than the angular frequency of the noise sources. 



1. Introduction 



^^ Turbulent eddies convecting within the boundary layer of an aerofoil are scattered 

^-^ into sound at the trailing edge. This, turbulent boundary layer trailing edge noise, 

QQ is a major source of broadband noise for an aerofoil in a uniform flow. It is the 

.^ dominant noise source for large wind turbines (Oerlemans and Sijtsma, 2007). For 

(^ installed fans, trailing edge noise corresponds to the minimum achievable noise 

cn level (Wright, 1976). 

^— I Although this paper focuses on trailing edge noise, rotors and propellers 

►>. are subject to several other noise sources (Hubbard, 1991). Some of these are 

•'"] broadband, for example leading edge noise (Paterson and Amiet, 1982; Homicz 

rN and George, 1974) due to upstream turbulence, tip vortex induced noise and stall 

^ noise (Moreau et al., 2009). Others are tonal, as in the case of: steady loading in 

the reference frame of the rotor (Gutin (1948), Hanson (1980)); periodic unsteady 

loading produced by stationary distortions in the flow. 

Trailing edge noise can be predicted in the time domain by solving the 

Ffowcs- Williams & Hawkings (FWH) equation (Ffowcs Williams and Hawkings 

* Author for correspondance (s.sinayoko@eiig.carn.ac.uk) 



Proc. Roy. Soc. A. 1-28; doi: 10. 1098/rspa. 00000000 
May 1, 2013 

This journal is © 2011 Tiic Royal Society 



(1969); Morfey and Wright (2007); Najafi-Yazdi et al. (2010)), as demonstrated 
by Casper and Farassat (2004). However, this type of prediction requires numerical 
differentiations and the calculation of retarded times that are computationally 
expensive. Most importantly, the roles of blade geometry and operating conditions 
become clearer with the frequency domain formulations that are the subject of 
this paper. Another key advantage of the frequency domain formulation is that it 
is more suited for a statistical description of the noise sources. 

One of the most successful frequency domain formulations was developed by 
Amiet (1976; 1978). Assuming knowledge of the pressure fluctuations travelling 
towards the trailing edge, Amiet 's formulation gives an analytical expression for 
the power spectral density in the far fleld. This analytical expression makes it 
efficient and attractive for fast turn-around applications. It has been applied 
to numerous applications including low speed fans (Rozenberg et al. (2010)), 
helicopters (Paterson and Amiet (1982); Amiet (1986)) and wind turbines (Glegg 
et al. (1987)). Although Amiet's approach was initially restricted to high 
frequencies relative to the chord {kC>l), it has been extended to lower 
frequencies {kC > 0.1) by Roger and Moreau (2005). 

The effect of rotation on trailing edge noise modelling was analysed in Amiet 
(1977) and Schlinker and Amiet (1981) for a rotor in a uniform flow. The flrst step 
is to estimate the instantaneous power spectral density radiating from the blade, 
while it is located at a particular azimuthal angle around the hub. Secondly, the 
power spectra are averaged in time. The final expression involves a Doppler factor 
of the form {oj' /oj)"', where uj' is the source frequency and oj the frequency of the 
observed sound. The exponent a takes the value 1 in Amiet (1976); Rozenberg 
et al. (2010), 2 in Schlinker and Amiet (1981) and -2 in Blandeau (2011). It is 
clear that an independent formulation is needed to explain these discrepancies 
and identify the correct form of Schlinker and Amiet's theory. 

This paper presents a new formulation of trailing edge noise for small rotating 
blade elements (section 2). It is analogous to the formulations developed by Glegg 
and Jochault (1998) for ducted fans, although cascade effects are not taken into 
account. The derivation starts with Goldstein's (1976) expression for loading noise 
and uses the Green's function of Garrick and Watkins (1954), which incorporates 
the effect of mean fiow, to derive an expression for pressure in the frequency 
domain. This expression is then used to derive the instantaneous PSD and 
the time averaged PSD. The new formulations generalize the compact source 
solution of Kim and George (1982), who used the expression of Ffowcs Williams 
and Hawkings (1969) for a rotating point load. They incorporate an accurate 
description of the blade rotation, hence of the effect of blade acceleration; this 
effect is not taken into account in Schlinker and Amiet's formulation where 
the blade is assumed to be in uniform rectilinear motion. The time averaged 
formulation can therefore serve to explore the validity of Schlinker and Amiet's 
approach. 

The derivation of the new formulation was inspired by Blandeau and Joseph 
(2011) who first compared Schlinker and Amiet's model with a solution of 
the FWH equation. Their formulation of Schlinker and Amiet's model differed 
significantly from that of Schlinker and Amiet (1981) but was not accompanied 
by an alternative derivation. Also, they neglected the effect of mean fiow and 
spanwise wavenumber. Both effects are taken into account in our new formulation. 



We revisit Schlinker and Amiet's theory and present simple derivations for 
some key steps that help to obtain the correct result (section 4). Of particular 
interest is the derivation showing how a theory derived for a stationary aerofoil can 
be generalized to an aerofoil in uniform motion. Finally, the results are validated 
by comparing the new formulation to Schlinker and Amiet's for multiple test cases, 
including a wind turbine, a cooling fan and an aircraft propeller (section 5). 



Nomenclature 

(X,Y,Z) Blade coordinate system (see fig. 4 and A. 10) 

( )d Subscript denoting division by Doppler shift 1 — Mz cos Og 

^m,n Kronecker delta: 1 if tti = n, otherwise 

r] Exponential decay rate of the spanwise coherence function 

Vo Gradient relative to Xo 

L(x,t) Unsteady lift force per unit area on a fiat-plate aerofoil 

Mx Mach number of the air relative to the blade at 0° angle of attack 

Xo Observer position 

Xp Present source position 

oj Angular frequency of the sound at the observer location 

oj' Angular frequency of the source 

^qqikx, ks) Surface pressure spectrum close to the trailing edge 

^ Source azimuthal angle relative to the trailing edge 

^{kx, ks, kc) Acoustically weighted lift 

(T Phase radius (equation (2.4)) 

( ) Unit vector in rotor plane in ( )-direction 

C chord 

Co Speed of sound 

k Acoustic wavenumber oj/cq 

kijkc Chordnormal (equation (2.14)) and chordwise (equation (2.14)) acoustic 
wavenumber respectively 

kr,kz Doppler shifted acoustic wavenumbers defined below equation (2.9). 

ks Spanwise acoustic wavenumber (equation (3.5) 

Kx,Kr Chordwise and spanwise hydrodynamic wavenumbers 



Kro Defined below equation (2.28) 

Is Spanwise correlation length 

Mx Cliordwise Macli number of the flow 

P(Kx,Kr) Wall pressure in the wavenumber domain 

Rpp{t) Autocorrelation function of the pressure at the observer location 

S span 

s Amplitude radius (equation (2.4)) 

S' {ijo) Pressure PSD produced by a stationary aerofoil immersed in a flow and 
measured by a stationary observer 

Spp{uj) Pressure PSD at the observer location 

Sqq{(jj) Surface pressure frequency PSD at the trailing edge. 

Uc Convection velocity of an eddy in the boundary layer near the trailing edge 

Uz Flow velocity in the minus z-direction (figure 1) 

X Source position along the chord from the trailing edge [X < 0) 

2. New formulation of trailing edge noise for rotating blades 
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(a) Pressure field in the time domain 

Consider the problem of trailing edge noise radiation in a uniform flow from 
a stationary rotor towards a stationary observer. Although the results of this 
section will be applicable to propeller blades or fans, the rotor is described as a 
wind-turbine whereby the blade is being rotated by the incoming flow (flgure 1). 
To derive the pressure fleld in the frequency domain we start by expressing the 
pressure fleld in the time domain. From Goldstein (1976), loading noise can be 
expressed as 

rT r 

p{xo,t) = — 'VG{x,T\xo,t) ■ L{x,T)dT,dT, (2.1) 

J-T Js 

where —T <t <T is the interval of time over which sound is emitted (T will be 
allowed to become inflnitely large), S is the blade planform and G a free fleld 
Green's function that satisfles the convected wave equation 

V2 - 1-^") G{x, t\xo, t) = -5{t - t)6{xo - x), (2.2) 

where Do/Dt = d/dt — Uzd/dzo (with f/^ > in the minus z-direction as shown 
in flgure 1). Compared to alternative methods, a convected Green's function 
simplifles the derivation by taking care of sound waves convection by the flow. 
Garrick and Watkins (1954) derived an expression for G that is analogous to the 
free fleld Green's function, 

G(x,rK,t) = ^(^-----/^°\ (2.3) 

47rs 

where x, Xo, t and r denote respectively the (present) source position, the observer 
position, the reception time and the emission time, and where the phase radius a 
and amplitude radius s are deflned as 

a={M,{z^-z) + s)/p\ s^ = {z^-zf+P\yo-yf + {zo-z)\ (2.4) 



where /3 = ^Jl^^M^. The geometric interpretation of a and s (Garrick and 
Watkins (1954)) is given in flgure 2, for the equivalent problem of a rotor and 
observer moving at Mach M^z in a quiescent medium. 

The Green's function takes a simple form in the far fleld (Hanson (1995)) when 
the source position is expressed in cylindrical coordinates (r, 7, z) and the observer 
position in spherical emission coordinates {Rg, 9^,^0)1 where Re corresponds to a 
in flgure 2 for x = (present source at the hub) : 

a ss i?e — -Zd cos^e — ^D sin6'eCos(7 — 7o), s^ Re{l — Mz cos 6e), (2.5) 

where the D subscript means that the variable has been divided by the Doppler 
factor 1 — MzCOsOe- The above expressions can be obtained by applying Taylor's 
approximation to equation (2.4) about the origin of the \x\ coordinate system, 
retaining terms up to order for s and 1 for a, and by converting the observer 
position to emission coordinates. 



(6) Pressure field in the frequency domain 

Since Amiet's trailing edge noise theory is expressed in the frequency domain, 
we seek a general expression for pressure in the frequency domain. Taking the 
Fourier transform over observer time t of (2.1) yields 



p{xo,oj) 



-T 



VG(a;, t\xo, uj) ■ L{x, r) dT, dr, 



(2.6) 



where the Fourier transform pair (/, /) is defined as 
1 



/(^) 



2tt 



me 



—iujt 



dt. 



fit) 



f{uj)e''''duj. (2.7) 



The frequency domain Green's function G is obtained by taking the Fourier 
transform over t of (2.3) and incorporating the far field approximation (2.5). The 
result is further decomposed into a series of azimuthal modes (see appendix 1) to 
simplify the computation of the gradient, 



+ 00 



Gix,T\xo,u;)=g{u;/co,Re,e,) J] J„(A:,r)e-^"(^o-'^/2)g-.(.r-n7-fc..)^ (2.8) 



g{k,Re 



'ikRa 



Sn'^Reil -M.CGSI 



(2.9) 



where we have introduced the axial wavenumber kz = ko cos 9e and the radial 
wavenumber kr = k^ sin Oe- 

From figure 3(b), the lift lies in the plane (7, z) and can be expressed in 
terms of the magnitude L and the pitch angle a. The dot product L ■ VG can 
be expressed in terms of the gradients of G in the azimuthal and axial directions 
that are easily derived from (2.8). Substituting the expression for the dot product 
into (2.6) yields a harmonic pressure of the form 



p{Xo,U}) 



W 



+00 



Le-^('^^-"T)drdS, (2.10) 



-T 



where /c/, = kz cos a — (n/r) sin a is the magnitude of a modal wavenumber 
orthogonal to the blade planform. 

We now turn to the rotational motion of the source. We will assume that 
the pitch change axis (PGA) of the blade is at the trailing edge (TE) and that 
the trailing edge alignment and face alignment (Hanson (1980)) are 0. In other 
words, lean and sweep are both neglected. The geometry of the blade is illustrated 
in figure 3. From figure 3(b) the azimuthal angle of a source within the blade 
planform is 7 = ilr + ^. If T tends to infinity, the time integral in (2.10) reduces to 
2ttL{x, ijJ — nQ)e'^"''^: the effect of the blade rotation is to select, for each azimuthal 
mode, the frequency uj — nil from the spectrum of the blade loading, so 



p{Xo,u}) = 2-Kig 



+00 

y^ g-m(7o-7r/2) 
i=— 00 



kLJnL{x,io - nQ)e'^''^'+''^Uj:. (2.11) 
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Figure 3. Blade geometry and source position. The blade path is unwrapped in the azimuthal 
direction -y: it moves at speed rfi along the -y-axis. The pitch change axis (PCA) is located at the 
trailing edge (TE) and indicates the axis around which the blade surface is rotated through pitch 
angle a relative to the rotor plane (r, 7). The source position is (r, f2r + i/), z) in the hub-fixed 
cylindrical coordinate system (figure 1) and (r, X) in the blade-fixed coordinate system. 



Finally, the lift L{x,uj) is given by Amiet's theory (detailed in section 3): 
for every spanwise wavenumber Kr, a gust of amplitude PQ{oj,Kr) travelling at 
convection speed Uc and scattered at the trailing edge, gives rise to a pressure 
jump PQ{uj,Kr)go{X,uj/Uc,Kr) along the chord, where go is the blade response 
function (Amiet 1976, 1978) and X is the chordwise source position measured 
from the trailing edge {—C<X<0 as shown in figure 3). Following Amiet, the 
gust amplitude Pq is a function of angular frequency oj and radial wavenumber Kr- 
Here, we express the amplitude as a function P{Kx, Kr) of chordwise wavenumber 
and radial wavenumber. Since the chordwise wavenumber of the gust equals oj/Uc, 
it is easy to show that P{u}/Uc,Kr) = UcPo{uJ,Kr)- The blade loading due to 
trailing edge noise is therefore given by 



L(x,uj) 



"+00 1 

—P{uj/U„ Kr)go{X, co/Uc, K.)e-^^'-'' 

-00 ^c 



dKr 



(2.12) 



As illustrated in figure 3 , the source coordinates {ip, z) can be expressed in 
terms of X as 

rip = —X cosa, z = —X sma. (2.13) 

The phase kzZ + nip = —kcX, where kc is a chordwise acoustic wavenumber: 



kc 



cos Q — sm a 
sin a cos a 



kz 
n/r 



(2.14) 



Note how the acoustic wavenumbers in blade coordinates {ki^kc) are obtained 
by rotating k^z + {n/r)^ by angle a around the PCA. After expressing the lift 



in terms of Amiet's blade response function using (2.12) and since dS = dXdr in 
(2.11), we find that the impact of the Uft on the pressure field is measured through 
the acoustically weighted lift ^ l defined as 

o rO 



^L{Kx,Kr,kc) 



C 



-c 



90 



{X,Kx,Kr)e-"'^''dX, 



(2.15) 



for which an analytical expression will be provided in equation (3.6). The 
frequency domain pressure at the observer location Xo can be expressed as 



p{Xo,(. 



i exp{—ikRe 



+ 00 

E 



-oo '-^c 



-J„(A;^r)e-^''(^o-^/2) 



SnReil -M^cosOe _ _ 

X PiKx,Kr)^LiKx,Kr,kc)e~'^-''dKrdr, (2.16) 



where Kx = (w — nQ)/Uc- 

(c) Instantaneous spectrum 

We now seek an expression for the instantaneous (Wigner-Ville) 
spectrum (Flandrin, 1998) for direct comparison with Schlinker and Amiet's 
theory (equation (4.11)). An autocorrelation function for the pressure p at the 
observer position Xq can be defined as 

Rppixo, t, t) = E[p{xo, t + Tlt)p\xo, t - t/2)], (2.17) 

where E denotes the expected value and • the complex conjugate. Since the blade 
element is rotating around the axis, this function is periodic in t with period Tn- 
The instantaneous spectrum Spp{Xo,uj,t) is defined as the Fourier transform of 
Rpp{xo, t, t) over r. An equivalent definition is 



bpp[Xo, ^, t) 



Rpp{Xo,io,zu)e™ dw, 



(2.18) 



where Rpp{Xo, oj, w) = E[p*{Xo, w — w/2)p(Xo, oo + tZ7/2)] is a cross-correlation in 
the spectral domain. Equation (2.18) states that the instantaneous spectrum is 
the inverse Fourier transform over frequency vj of Rpp(xo,oj,iu). 

Spectral cross-correlation Rpp{Xo,UJ,w) 

From (2.16), the spectral cross-correaltion is given by 



RppiXo,OJ,w) 



[8vri?e(l-M,cos^e)]' „i-^ J^^ 



gi(n-n')(7o-'r/2)^ 



■Rt 

Rh 



■Rt 

Rh 



+ 00 



E[P*{Kx, Kr)P{K'x, K)] dKl. dKr dr dr, (2.19) 



where the primed wavenumbers and non-primed wavenumbers are associated 
respectively with frequency a; + tZ7/2 and uj — w/2. 

We will now derive two approximate expressions Rpp and Rpp for Rpp, by 
simplifying the integrals over r' and r respectively. These two expressions will 
then be combined to give an expression for Rpp that has Hermitian symmetry; 
this is necessary for the instantaneous PSD to be real. 

From Blandeau (2011) (appendix D), if we assume that the flow is 
approximately 2D and that it is uniform across the span (strip theory assumption) 
then the two gusts are correlated only when their wavenumbers are equal: 
Kr = K'^ and Kx = K'^ (i.e., from the definition of Kx below equation (2.16), 
\i w = {n' — n)Q), and we have 

E[P*{Kx, Kr)P{K'x, K)] = U', 6{K', - Kr) 5 {w - {n' - n)Q) %g{K'x, K), 

(2.20) 
where ^qq is the wavenumber spectrum of the aerofoil surface pressure produced 
by the turbulence. The above equation assumes that the convection velocity is 
constant over the area where P{Kx, kr) and P{K'^, k'^) are correlated. For a given 
radius r, that region is defined by r — Ar <r' <r + Ar, where Ar is slightly larger 
than the turbulence correlation length Is- Since f7^ = r'fi, variations in convection 
velocity are negligible over that region provided that Ar <C r', i.e. Is ^ r' . Using 
Corco's equation (3.4), at high frequency, 1$ ~ 0.8r'r2/(aj'?7) (where r/ is a constant 
of order one) so the necessary condition becomes i7 ^ w': the rotation speed should 
be much smaller than the source frequency. This corresponds to the assumption 
made by Amiet (see introduction of section 4). 

Substituting (2.20) in (2.19) leads to Rpp = Rpp, where 

fRt +°° +°° 

•'^h n=~oo n'=—oo 



5{w- {n - n)Q) dr, (2.21) 



B{w)= ^^ ^^^, (2.22) 



g-jro_Re/co 



'+00 

I',{u;,ru,r)= Ln,n'{K)^ggiK'x,K',) dK',, (2.23) 



Ln,n'iK'j.) 



■^* kLCk'^C 



n^^'^Jn (krr) Jn' {Kr') e^^^(---') dr'. (2.24) 



The integration over r' may be truncated at r' = r it Ar. In this narrow strip, 
the flow variables can be considered constant in the amplitude terms of (2.24) 
(fci, /c^, C/c) t^c) C' and C): they are equal to their values at r' = r. We seek to 
simplify (2.24) and put it in an exponential form to capture its phase. This is 
needed to analyse the spanwise behaviour of the solution. The main issue is with 
^L, which may be studied using a power series expansion. However, for simplicity, 
we will neglect the phase changes in ^x, and assume ^^(r') ~ ^'x^ir) so that L^n' 
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takes the form 



Ln,n'\Kj., 



) n'i 



~~Ur. 



■niKx,K,kc)^LiK'x,K',,k'c)Jn {Kr) Pn,n'iK), 



where, making the change of variable ri = r' — r, 



J~^n,n' (,-^r^ 



rAr 



-Ar 



Jn' {k'Av + r)) e-'^'rV dr^_ 



(2.25) 



(2.26) 



The Bessel function in (2.26) can be rewritten as an integral (Abramowitz and 
Stegun (1972)): 



Jn'iKiv + r)) 



1 

2^ 



^i{n''y~k'^{r]+r) sin7) 



dj. 



Substituting (2.27) into (2.26), 



J^^n,n' \^r/ 



1 

2^ 



gj(n'7— fc^r sin7) 



+Ar 



-Ar 



<K~Koyndr,dj, 



(2.27) 



(2.28) 



where K'j.q = —k'^ sin 7. If Ar is large relative to the wavelengths of interest in this 
problem, i.e. if K'j.Ar S> 1, then the second integral in (2.28) simplifies to 

f+00 

^-iiK'^-K'^^)^ dr? = 27r6{K', - K',q), (2.29) 

DO 

The above simplification may not always be possible in which case the integral 
can be solved numerically. Substituting (2.29) into (2.28) to get PnyiK^.), and 
combining the result with (2.25) and (2.23) yields 



iLn'i^^'^^r) 



-^ Jn [Kr] Qn'^ni^, ^, r) 



(2.30) 



where 



J(n'j-k'^r sin 7) ^,* 



n(Kx,Ko, kc)^LiK'x,Ko, k'c)%,{K'^, K',^) dj. 



(2.31) 



?(2) 



We can derive a different expression Rpp for Rpp, by approximating the integral 
over r in equation (2.19) and by following the steps in equations (2.21) to (2.31) 
to get 



R^I^{xo,uj,w) = B{w) 



^h n=~oo n'=~oo 



5{w- {n - n)n) dr, (2.32) 
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where 

kjk' C^ 
/„/,„(a;,ro,r')~ — -^ — J„, (/c^r') Q„^„/(w, tz7,r'), (2.33) 

J— TT 

(2.34) 

where KrQ = —krSm.^ and the wavenumbers fc^, /c^, kc and /c^ are evaluated 
at radius r'. Note that in deriving K^o we have used the conjugate form of 
equation (2.27) to express Jn{krr). We have also chosen this time Kx instead 
of K'j^ as the first argument of ^qq. This is justified since Kx equals K'j^. 

If we change ro to — ro and exchange n and n' then Q'^, ^ turns into Qj^ „/. This 

is sufficient to show that Rpp {Xo,uj, —w,r) equals the conjugate of Rpp {xo,io,'oj). 
We express Rpp as 

^PP^ ""t "' . (2-35) 



which satisfies the Hermitian symmetry required for the instantaneous spectrum 
to be real valued. 

Instantaneous spectrum Spp{xo,oj,t) 

Substituting equations (2.21) and (2.32) into (2.35) and then into (2.18), 



bpp(Xo, U), t) 



[■Rt +2f 

B{mQ)e 

Rh 



V R('"7.nv-*™(')'°-'^/2)- 



m=—oo 



Y. /„,„+„(a;,mJ^,r) e^^^dr, (2.36) 

\n=— oo / 

where we have made the change of variable m = n' — n and 

Note that the expression of equation (2.36) is always real but it can become 
negative due to a known defect of the Wigner-Ville spectrum, which does not 
define a true spectral density. 

Time averaged PSD Spp{Xo,uj) 

The time averaged PSD is obtained by retaining the mode tti = in (2.36), so 
w = and the prime and non-prime quantities are at the same frequency w, and 
since /* „ = /^ „, we can write 



jpp^Xo, (jJ) 



Rt +°° 



Rh 



B{0) y iniuj,r)dr, (2.37) 



n=~OD 
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where, from (2.30) and (2.33), 






COS (n7 — krr sin 7) x 
|^(i^x, -kr sin7, %)|2$gg(Kx, -A;^sin7) ^7. (2.38) 



Our numerical experiments using the above equations indicate that Spp is always 
positive. 

For a small blade element of span S, 

+00 
Spp{xo,oj) = B{0)S Y, ^n(^,r). (2.39) 



3. Trailing edge noise theory for isolated aerofoils 



Observer 




Figure 4. Coordinate system for an isolated aerofoil in a wind tunnel. The aerofoil and observer 
are stationary. The flow is uniform at zero angle of attack and Mach Mx ■ 



Consider a flat plate in a uniform flow of Mach number Mx at zero angle of 
attack (flgure 4). The observer location is expressed using a cartesian coordinate 
system {X,Y,Z), with X in the chord wise direction and pointing downstream 
and Z in the vertical direction. From Amiet (1976, 1978), using the large span 
assumption (equation (2.29)), the PSD at frequency lo is given by 



Spp(U}) 



uj C Z 



CO 2 27rs2 J 2U, 



S vr 



'^qq(.kx,ks)\'ifL(.kx,ks,kc)(' 



(3.1) 



where 



x2 + /?2(y2_^z2). 



/?2 = 1-Ml, 



Uc = 0.8MxCo, (3.2) 

and from Roger and Moreau (2005) the wavenumber spectrum and spanwise 
correlation length are given by 

vr 



U, 



^qq{kx, ks) = ls{kx,ks)Sqq{u}), 



(3.3) 
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The term rj represents the exponential decay rate of the spanwise coherence 
function; Brooks (1981) measured r/ = 0.62 for a NACA 0012 at Mach 0.11 and 
zero angle of attack. For simplicity, this value will be used for all Mach numbers 
in section 5. For a stationary fiat plate in a uniform flow (Amiet (1975)) 

kx = ^, ks = --, kc = ^(Mx--). (3.5) 

Uc Co s co/3^ V s J 

(a) Acoustic lift 
The acoustic lift, defined in equation (2.15), can be expressed as 

^L{kx, ks, kc) = jl y^.^ eii[^2{iB - iA)) + e^^^ [l - ei{{V2iB 

(3.6) 

A = kx + kc, B = kx + K + MxJl, fi = Mckx/P^, (3.7) 

where the overbar denotes normalisation by C/2 and Mc = Uc/cq the Mach 
number of the turbulent eddies in the boundary layer close to the trailing edge. 
The square roots in equation (3.6) have a branch cut along the negative real axis. 
The wavenumber k is a function of the spanwise wavenumber ks- It is defined 
as 

/x Jl - fc|/(/3^)2 if kl<{Piif 

(3.8) 



^-i\^\^^/{j^i^l if kl>{pfi) 



2 



SO that the imaginary part of n is always negative. This is required for the error 
functions in (3.6) to converge in the far field. The square roots in equations (3.8) 
and (3.6) are classicaly defined with a branch cut along the negative real axis. 

Figure 5, shows how |^| varies with spanwise wavenumber, for a constant 
frequency of 1 kHz (solid line). The transition between supercritical and 
subcritical gusts occurs at ks/ P^i = 1. It can be seen that a sharp increase occurs 
near ks/l3iJ.= 1. That increase is non-physical and occurs because the governing 
equation reduces from a Helmholtz equation to a Laplace equation. This was first 
pointed out by Roger and Moreau (2005). Since ^ is asymptotically constant 
at low and high values of ks we define ^ as a piecewise function, such that it 
is constant for supercritical gusts {ks//3iJ.< 1), and for highly subcritical gusts 
{ks//3fi > 10). Between those values, we use a linear interpolation of |^| in terms 
of logiQks/l3fi. The piecewise implementation of |^p is shown as a dashed line 
in the figure 5. 

(6) Surface pressure power spectral density 

The PSD of the incoming pressure fluctuations at the trailing edge, denoted 
Sqq{uj), can be measured experimentally. If no experimental data is available, it 
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Figure 5. Acoustic lift spectrum \^{kx, ks, fec)|^ as a function of spanwise wavenumber ks, at 
Mach Mx =0.1 and frequency / = 1000 Hz. The radiating wavenumber kc is defined assuming 
an observer at x/s — 0.7. The transition between supercritical and subcritical gusts occurs at ^ = 
ks/iPfJ,) ~ 1. The acoustic Uft based on equation (3.6) becomes unphysical near the supercritical- 
subcritical transition (solid line, see Roger and Moreau (2005)), so we assume that |'I'| is constant 
for ^ < 1 and ^ > 10, and use a linear interpolation as a function of logj^Q (5) for 1 < C ^ 10 (dashed 
line) . 



can be estimated by using empirical low-order models. These low order models are 
expressed in terms of parameters characterising the boundary layer, such as the 
boundary layer thickness, displacement thickness, wall shear stress etc. A review 
of the different models is given by Blandeau (2011). This paper uses the model 
of Chou and George (1984) which gives 

^.<?M = (^pf^i) ^^(^)' (3-9) 

where S* is the boundary layer displacement thickness, Ux the chordwise flow 
velocity and uj = uj6* /Ux, 

6* 



rC7(24.3 + 0.6625x)10"^ if x < 4° 

C(26.95 + 0.6625(x - 4) + 0.3044(x - 4)^ + 0.0104(x - 4)^)10"^, x > 4°, 

(3.10) 



1 



where x is the angle of attack and 

1.732 X IQ-^cJ 
1 - 5.489W + 36.74a;2 + 0.1505^5 

1.4216 X lO^-'^tJ 



F(lJ) 



if ZU<0.06, 



L 0.3261 + 4.1837aJ + 22.818a;^ + 0.0013aJ'^ + 0.0028a; 



?7tt3 



i7:^5 



if aJ>0.06. 
(3.11) 



(c) Discussion 

A number of approximations were made in this section to simplify the 
comparison between the formulation of section 2 and that of Amiet (section 4). 
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These approximations may need to be replaced with more accurate models for 
practical applications. 

The first approximation is the equation (3.4), giving the correlation length Is, 
that was derived by Corcos (1964) for a turbulent boundary layer over a flat plate 
with zero pressure gradient. A more general equation has been proposed by Roger 
and Moreau (2005). 

A second approximation is the spanwise dependence of the acoustic lift 
described in flgure 5. This approximation was made for convenience and does 
not affect the rest of the analysis. 

Thirdly, the wall-pressure spectrum of Chou and George (1984) in section 3(b) 
may be replaced by more complex models such as that of Rozenberg et al. (2010). 
Rozenberg's model takes into account the influence of an adverse pressure gradient 
and should therefore better predict the effect of aerofoil camber and thickness. 
However, more research is needed in this area and, ideally, the wall-pressure 
spectrum should be measured experimentally. 



4. Schlinker and Amiet's approach for rotating blades 

This section is a detailed and corrected derivation of the theory presented 
by Schlinker and Amiet (1981). In particular, the derivation of the present source 
position is both simpler and more general. This section explains how to apply a 
theory derived in the wind tunnel reference frame, where a stationary source and 
observer are immersed in a uniform flow (see section 3), to the general case, where 
both the source and the fluid are moving relative to the observer. 

Amiet (1986) observed that the sound emitted by a rotating blade is 
approximately equal to that emitted by a translating blade at the same location. 
His observation was based on the expression of Lowson (1965) for the pressure 
radiating from a rotating dipole: the difference between a rotating blade and a 
translating blade lies in the acceleration of the source in the direction of the 
observer. Amiet argued that when the angular velocity 0, is much larger than the 
source frequency oj', the acceleration term becomes negligible. This high frequency 
assumption will be validated in section 5 against the new formulation of section 2. 
In the following, we will assume that the blade is in rectilinear motion. 

In the situation we are considering here, the observer is stationary relative to 
the hub and the blade moves through the fluid (flgure 1). However the formulae 
of section 3, for a blade in a uniform flow, assume that the observer is stationary 
relative to the blade. We therefore examine the current problem successively in 
these two reference frames. 

In the reference frame of the observer (O), the observer is stationary and the 
blade moves rectilinearly at Mach Mbo in a- uniform flow of Mach Mpo- Consider 
a pulse of sound emitted at Xe from the blade. When this pulse is received at the 
observer position Xo, the blade has moved to the present source position Xp, as 
illustrated in flgure 6(a). 

In the reference frame of the blade {B), the blade is stationary and the source 
is flxed to the present source position {x^ = x^). The Mach number of the flow is 
Mfb = Mpo ~ Mbo a-nd the observer is moving at Mach Mqb (flgure 6(b)). 
Comparing flgure 6(b) with flgure 4, this situation is equivalent to that of an 
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MboCqZ 



MpoCaT, 



(a) Reference frame of the observer (O): stationary observer, blade moving in uniform fiow. 






MfbCoT, 



MobCqT,, 



(b) Reference frame of the blade (B): equivalent to isolated aerofoil case of figure 4 if Mob ~ 0. 




MBFCoTe xi 

(c) Reference frame of the fluid (F): stationary fluid so \x^ 



'CoTe 



Figure 6. Noise radiation from a blade in uniform rectilinear motion analysed in three different 
reference frames that are fixed to: (a) the observer (O); (b) the blade (B); (c)the fluid (F). 
Mab stands for "Mach number of A relative to B", where A and B correspond to one of the 
reference frames. All figures use the same scale so Mach numbers are comparable. In each case, 
a pulse of sound is emitted from the emission position Xe at emission time and reaches the 
observer Xo at reception time, while the source moves with the blade to the present position Xp. 
The convected source position Xc is obtained by convecting the emission position at the Mach 
number of the fluid for the duration of propagation Te- The triangle between present source 
position, convected source position and observer position, shown in dotted lines, is independent 
of the reference frame: the vectors Xo — Xp, Xo — Xc and Xp — Xc are invariant. The position 
of the blade is shown in grey at emission time and reception time, assuming for simplicity that 
MpB is aligned with the blade chord (zero angle of attack) . 
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isolated aerofoil in a uniform flow provided that the observer is stationary, i.e. 
Mob = 0. Here, the observer is moving relative to the blade which leads to a 
Doppler shift whose impact on the instantaneous PSD will be discussed later. 

The above analysis indicates the steps required to apply the formulae derived 
for an aerofoil in a uniform flow to the present problem: 

• move the origin of the coordinate system from Xe to the present source 
position Xp] 

• rotate the axes of the coordinate system if necessary; 

• account for the movement of the observer relative to the blade (Doppler 
shift). 

Present source position 
From flgure 6(a), 

Xp = Xe + MsoCoTe, (4.1) 

where Tg denotes the propagation time. The above result is more general than the 
equation (46) in Schlinker and Amiet (1981) that is valid only if the chordwise 
Mach number of the flow equals MpB, i-e. if the angle of attack equals zero. 
The source position at emission time satisfles \xe\ < Rt, where Rt is the blade tip 
radius, while Tg ^ 1 for propagation to the far-fleld, therefore (provided M^o 7^ 

0) 

Xp^MBoCoTe- (far held). (4.2) 

For an observer in the far field, we can hence assume a^e = so that the source is 
located at the hub. 

Propagation time 

As demonstrated in the following, the propagation time Tg is related to 
the convected source position Xc- The convected source position is defined as 
the emission position x^ in the reference frame of the fiuid (figure 6(c)). In 
that reference frame, the fiuid is stationary so the wavefronts are spherical and 
acoustic waves are propagating at the speed of sound. This allows us to relate the 
propagation time Tg to Xo and Xc since 

where we have used the fact that Xo — x^ is independent of the reference frame 
(see figure 6). 

In general the convected source position Xc is obtained by seeding the fiow 
with a small particle that is released from the source position at emission time 
and convects at the Mach number of the ambient fiuid; Xc corresponds to the 
position of the particle at reception time. Using the reference frame of the observer 
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(figure 6(a)), we get 

x^ = x^ + MrocoTe « MpoCoTe, (far field). (4.4) 

Substituting (4.4) into (4.3) and taking tlie square of tlie result, 

{coTef^\xo-MFoCoTe\^, (far field). (4.5) 

The above equation leads to a second order polynomial equation in R^ = c^Te 
whose solution is given by 



R ( -Mfo cos e + y'I - Mj^c, sin^ 9 
Re = -^ ^ —^ ' , (far field) , (4.6) 

where denotes the angle between Mpo and Xo- 

Coordinate system rotation 



blade movement 




\ 

wind \ a 

\^ blade movement 


\ 


^ 


X 


A^ 






wind 

(8) 



(a) (xi,y2,Z2) is obtained by rotating the (b) {X,Y,Z) is obtained by rotating 
{xi,yi, zi) by n/2 — 7 around the zi-axis. (x2,y2, Z2) by q around the 2y2-axis. 

Figure 7. Coordinate systems {xi, yi, zi), centred on the present source position Xp, {x2,y2, Z2), 
such that 2/2 is in the spanwise direction, and {X, Y, Z) , such that Y is in the spanwise direction 
and X in the chordwise direction. 

In the reference frame of the blade used in section 3, the y-axis is pointing 
spanwise and the x-axis is pointing chordwise. Let x-j_ = {xi,yi,zi) denote the 
coordinate system centred on the present source position {xi =Xo — Xp). For the 
y-axis to point in the spanwise direction, the a^i-coordinate system is rotated by 
7r/2 — 7 around the zi-axis (figure 7(a)), i.e. 



X2 = Rz{ti/2--i)x-s_, 



(4.7) 



where Rz denotes the rotation matrix about the zi-axis (see equation (A. 3)). 

Similarly, for the x-axis to point in the chordwise direction, the a;2-coordinate 
system is rotated by a around the 1/2-axis (figure 7(b)): 

X = Ry{a)x2, (4.8) 



where Ry denotes the rotation matrix about the y2-syi\s (see equation (A. 3)). 
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The transform from the observer reference frame (figure 6(a)) to the blade 
reference frame is hence given by 

X = Ry{a)RJy'K/2-a){xo-Xp). (4.9) 

Doppler shift 

The sound frequency cv at the observer is shifted compared to the emission 
frequency co' . The ratio oj/uj' is called the Doppler shift and is a function of 
the source Mach number Mbo smd the flow Mach number Mpo- As explained 
by Amiet (1974), the Doppler shift is given by 

^ = 1 + MboCO ^ (for field), (4.10) 

<^' 1 + {Mfo - Mbo) • CO 

where CO is the unit vector from the convected source position {X(. in figure 6) 
to the observer position. 

It can be shown (Amiet (1974) and appendix 1) that in the far field the PSD 
for a moving observer, Spp, is related to the PSD for a fixed observer, S' by 
Spp{u}) = {oj' /uj)S'pp{uj'), therefore 

Spp{xo,t,uj) = —S'{X,T,uj'), (4.11) 

(jj 

where S' can be estimated using the isolated aerofoil expression of (3.1). 

Furthermore, equation (3.1) is expressed in terms of the chordwise Mach number 

Mx = MpB • X of the flow relative to the blade. A derivation of Mx, extending 

that of Schlinker and Amiet (1981) to the case of a non-zero angle of attack, is 

given in appendix 1. 

The time averaged PSD is obtained by averaging the instantaneous PSD, i.e. 

equation (4.11), over one rotation of the rotor. Since the instantaneous PSD is a 

function of source time r, the time increment dt must be expressed in terms of dr. 

If n is the number of acoustic periods measured at the observer location during 

dt then dt = nlix juj. The time taken for the source to generate those n periods 

is dT = Tilix I Lij' so dt = [uj' /oj)dT. Finally, it is convenient to express dr in terms 

of the blade azimuthal angle 7 using the relation Vt = dj/dr. The time averaged 

PSD from a rotating blade element is hence given by 



^ r27r /, ,/\ 2 

zvr 



) S;piX,u;',^)dj. (4.12) 



The above result agrees with Schlinker and Amiet (1981), provided that oj is 
replaced by lOq in the left hand side of (54) in their report, and lo and lOq are 
swapped in the left hand side of their equation (56). Note that there has been some 
discrepancy in the literature about the exponent in the Doppler term of (4.12). For 
example, Amiet (1977) initially proposed a value of 1 for the exponent. The same 
is considered by Rozenberg et al. (2010). This ignores the additional weighting 
of the time increment and is equivalent to averaging over the angular position 
of a blade segment. We propose that the correct value of the exponent is 2. In 
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section 5, evidence for this is provided by comparison of predictions using the two 
methods presented in sections 2 to 4. 



5. Results 

In this section, both Amiet's formulation (equation (4.12)) and our new 
formulation (equation (2.39)) are applied to several model blade elements. These 
blade elements have been chosen to cover the range of applications featured in 
the literature on trailing edge noise: cooling fans (Rozenberg et al., 2010), open 
propellers (Pagano et al., 2010) and wind turbines (Glegg et al., 1987; Oerlemans 
and Schepers, 2010). Each blade element is described through its distance from the 
hub (radius), its geometry (chord and pitch angle a), and the set of Mach numbers 
(blade Mach MbOi Ao'^ Mach numbers relative to the observer Mz = Mfo a-nd 
to the blade Mx = Mfb), whose values have been chosen to match a typical 
situation. For example, a typical wind turbine element has a large chord, low 
pitch angle and rotates at low Mach number in a low wind speed. The specific 
values chosen here, presented in table 1, were the ones used by Blandeau and 
Joseph (2011). The span is arbitrarily equal to the chord. Note that for each blade 
element, the values are not independent; the three Mach number form a velocity 
triangle wherein the pitch angle is between Mpo a-nd Mfb- Furthermore, the 
values would usually change along a full blade. 

Cooling fan Wind turbine Open propeller 



radius 


0.40 m 


29 m 


1.80 m 


chord 


0.13 m 


2 m 


0.31 m 


Mbo 


0.0525 


0.165 


0.748 


a 


34deg 


10 deg 


13 deg (take-off), 38 deg (cruise) 


Mfo 


0.0354 


0.029 


0.584 (take-off), 0.228 (cruise) 


Mfb 


0.0633 


0.167 


0.949 (take-off), 0.782 (cruise) 



Table 1. Typical parameters for the blade segments used in three different applications: open- 
propeller, cooling fan and wind turbine. These parameters are the ones proposed by Blandeau 
and Joseph (2011). The span is defined as a third of the radius. The wind Mach numbers, relative 
to the observer Mfo{= M^) or to the blade Mfb{~ Mx) can be obtained from the pitch angle 
a and blade speed Mbo but are given for completeness. 



Although this paper focuses on noise radiation from small blade elements, 
the noise spectrum of a full blade may be obtained using strip theory. In strip 
theory, one divides the blade into small elements along its radius and sums 
the (uncorrelated) noise spectra generated from each element. Limitations and 
extensions of strip theory are discussed in Christophe et al. (2009). In addition to 
the approximations discussed in section 3(c), the thin-blade approximation has 
been used. This is not a good approximation for angles close to the propeller axis 
where a numerical approach taking the thickness of the blade into account may 
be preferable (Zhou and Joseph, 2006). 

Figure 8 compares the sound pressure level directivity obtained using the 
new formulation and Sclinker and Amiet's approach, when the exponents of the 
Doppler shift (see equation (4.12)) takes the values 1, 2 and -2 taken from the 
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(b) Open propeller at takeoff. 



Figure 8. Sound pressure level directivity using the new formulation (equation (2.39), dashed 
line) and Amiet's formulation (equation (4.12)) with the exponents of the Doppler shift set to 
—2 (o symbols), 1 (D symbols) and 2 (solid line). The frequency is fcC = 5. The decibel scale 
assumes jXo| = Im and wind speed is from right to left. 



literature. A low Mach number case (wind turbine) is shown in figure 8(a), and 
a high Mach number case (open propeller) in figure 8(b). Differences of up to 5 
dB and 20 dB can be seen when Schlinker and Amiet's approach are used with 
an exponent of 1 (square symbols) and —2 (round symbols) respectively. The 
best agreement between the new approach and Schlinker and Amiet's approach 
is obtained for an exponent of 2 (+ symbols), which validates equation (4.12). 

Figure 9 shows the sound pressure level directivity for the blade segments 
defined in table 1 that model a cooling fan (figure 9(a)), a wind turbine (figure 
9(b)), and an open propeller at take-off (figure 9(c)) and cruise (figure 9(d)). In 
each case, the results are obtained using the new formulation (dashed line) and 
Amiet's approach (solid line). In figures 8 and 9 the levels are normalised to give 
an effective observer distance of |a;o| = Im and the observer is located at 70 = 0°. 

The two approaches agree to within 1 dB in all cases. These results suggest that 
Schlinker and Amiet's theory is valid when the acceleration effects are negligible, 
i.e. cj ^> O, and when Amiet's blade response function is applicable, i.e. kC > 1. 
Since uj/Q, = kr/MsOi this ratio is larger than 2.9 for all the test cases in table 1. 
Note also that Schlinker and Amiet's formulation is applicable even for high speed 
open rotors: figure 9(d) shows a good agreement for the open propeller model at 
cruise for which the chordwise Mach number is very high {Mps = 0.95). 



Conclusions 



Schlinker and Amiet's (1981) theory of trailing noise from rotating blade has 
been validated by comparison with a new formulation that implicitly includes 
acceleration effects. The results presented in this paper show that the modified 
form of Schlinker and Amiet's predicts sound pressure levels within 1 dB of our 
new formulation for subsonic chordwise Mach numbers of up to 0.95 at high 
frequencies (fcC > 1 and oj^Vi). It is therefore applicable to both low speed 
applications (cooling fans, wind turbines) and high speed ones (propellers). 

We propose that the correct exponent for the Doppler weighting (oj' /co)"' 
in Schlinker and Amiet's theory is a = 2. Using a = l, as in Amiet (1977) 
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Figure 9. Sound pressure level directivity in decibels for the model blade elements of table f , 
and an observer Xo such that \xo\~lm (far field is implied) and 70 = 0°. The normalized 
frequency kC equals 0.5 (left column), 5 (middle) or 50 (right). Solid line is Amiet's approach 
(equation (4.12)), dashed line is the new formulation of equation (2.39). Wind speed is from 
right to left. 
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and Rozenberg et al. (2010) gives results that agree with our formulation to within 
5 dB. However, all other things being equal, putting a = — 2 in Schlinker and 
Amiet's theory overestimates the result by up to 20 dB at high Mach numbers. 

The range of validity of Amiet's theory for rotating blades is currently limited 
to high frequencies for two reasons: it uses a high frequency blade response 
function and neglects the effect of blade acceleration. More research is needed to 
identify the most important effect between the two. If the effect of acceleration can 
be neglected, i.e. at high frequency relative to the rotational speed, Amiet's theory 
can be extended to lower frequencies by using a low frequency blade response 
function (Roger and Moreau (2005)). 
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Appendix 

Far field Green's function 

Taking the Fourier transform of (2.4) over t, and using the far field approximation 
of (2.4) yields 

„-iui{T+s/co) p~iuJT 
G{x,t\Xo,Oj) = = S^^^+Kvcosi,^-^,,)) ^ ^p^-^^ 

where k^ = k^ cos 9^ and kr = k£, sin 9(,. Using the Jacob i- Anger expansion 

n=+oo 

gifc,rcos{7-7o)= Y^ J„(fc^r)e^"(^~''"+^/2)^ (A.2) 

n=— oo 

which yields equation (2.8). 



Rotation matrices 

(cos 9 — sin 6 0\ /cos ^ — sin 9\ 

sine cos6' , Ry{9)=[ 1 . (A.3) 

1/ \sin6l cos6l / 



Instantaneous frequency and Doppler shift 

Physical interpretation relative to the speed of time 

The Doppler shift indicates how the frequency of a pure tone varies for 
a moving observer (relative to the source), compared to a fixed observer. 
Following Amiet (1974), an alternative physical interpretation can be obtained 
by observing that, for a pure tone, the pressure p' for a fixed observer and p for 
a moving observer are given by 

f r,'(t) = e~*'^'* 

<^ so p(t) = e— '(-/-')* =p'((c^V)0' (A-4) 

[ pit) = e-*'^* 

where uj' is the source frequency and u the frequency received by the moving 
observer. The above equation implies that the measurement of the pressure time 
history is sped up (or slowed down), by a factor w/oj', for a moving observer, 
compared to a fixed observer. 
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Power spectral density in two reference frames 

From equation (A. 4), the autocorrelations Rpp{t) and Rpp{t) for a fixed 
observer and a moving observer (relative to the source), are such that 



Rpp{t) = lim 



1 



T-^+oo 2T 



1 



■T 



1 



p{T)p{t-T)dT= \\m 



■T 



-T 



p {aT)p {a{t — t)) dr 



raT 



lim ^ 

T-5>+oo 2aT 



p{T')p{at — t) dr' = R'{at) 



-aT 



(A.5) 
(A.6) 



where a = uj/oj'. 

Since the power spectral density (PSD) is defined as the Fourier transform of 
the autocorrelation, the PSD Spp{u}), for a moving observer, is related to the PSD 
Spp{uj') for a fixed observer by 



Spp{^) 



r+oo 



Rpp{t)e''^^dt: 



r+oo 1 

R!pp{at)e''^'''^ dt = -S'pp{J), (A.7) 

-^ Oj 



I.e. 



w 



Spp{u]) = —S'Juj'). 
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(A.8) 
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(a) Flat plate (in grey) in the (7, £) plane, 
where z is pointing against the flow and is 
orthogonal to the rotor plane, and 7 is in 
the direction of rotation. The angle between 
the rotor plane and flat plate is a and the 
angle of attack is X-The flow Mach number 
relative to the blade is Mfb = —M^z — 
M^-y. The chordwise Mach number is Mx- 



(b) Relationship between the flow Mach 
number M^ along 7, the blade Mach 
number Mbo and the (wind) cross-flow 
Me in the rotor plane {x,y). Note that 
the flat plate is not in the rotor plane 
(see figure (a)). The angle 7 is the 
azimuthal angle of the blade. The angle 6 
gives the direction of the cross flow relative 
to y. 



Figure A. 10. Coordinate systems and Mach numbers around a flat plate moving in a uniform 

flow. 



We derive the chordwise Mach number Mx given the blade Mach number 
Mbo = Mbo'1 ^^'^ the flow Mach number Mpo- As illustrated in figure A. 10(b), 
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Mpo can be decomposed into a component M^ normal to the rotor plane (inflow) , 
and a component Mq parallel to the rotor plane (cross-flow): 

Mfo = -M-^z - MeO. (A. 9) 

From figure A. 10(a), the chordwise Mach number Mx is given by 



Mx = MpB cos X = Jm^ + M| cos x- (A. 10) 

From figure A. 10(b), the azimuthal Mach number M^ of the flow relative to the 
blade can be expressed in terms of Mbo and Mq by 

M^ = Mbo + Mgcos{-f + e). (A. 11) 

Substituting equation (A. 11) into equation (A. 10), 

Mx = cosxVIMbo + Mg cos(7 + ^)]2 + Ml (A. 12) 



